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ELECTROMAGNETIC  SCATTERING 

This  interim  scientific  report  summarizes  the  work  carried  out  under  the 
Air  Force  Office  of  Scientific  Research  Grant  No.  77-3188(A)  during  the  year 
ended  31  December  1977. 

It  has  been  a year  of  steady  progress  on  several  fronts. ^One  of  the  major 
efforts  has  been  in  connection  with  the  design  of  broadband  absorbers  to  re- 
duce the  backscattering  cross  section  of  wing-like  structures  when  illuminated 
at  close  to  edge-on  incidence  with  magnetic  vector  parallel  to  the  edge.^^ 
Under  a previous  Air  Force  Contract  it  was  shown  that  thin  coatings  could  be 
simulated  using  a surface  impedance,  and  that  the  materials  available  could 
produce  a significant  cross  section  reduction  over  broad  aspect  and  frequency 
ranges,  but  the  time  available  did  not  permit  the  detailed  analytical  and 
numerical  investigation  necessary  to  determine  the  maximum  cross  section  re- 
duction attainable. 

Such  an  investigation  has  now  been  performed  under  the  present  Grant. 

Using  a 30®  angle  ogival  cylinder  with  a maximum  allowed  value  for  the  magni- 
tude of  the  surface  impedance,  analyses  and  computations  have  been  performed 
to  find  the  surface  impedance  profiles  that  are  most  effective  in  reducing  the 
backscattering  cross  section  for  H polarization,  with  particular  reference  to 
an  aspect  range  i 30°  about  edge-on.  Analytical  expressions  for  the  fields 
diffracted  by  the  edge  of  a uniform  impedance  wedge  and  by  a discontinuity  in 
the  value  or  first  derivative  of  an  impedance  on  a plane  surface  were  used  to 
specify  the  (local)  surface  impedance  to  minimize  the  direct  scattering  from 
any  singularity,  as  at  the  front  and  rear  edges  of  the  ogival  cylinder.  Various 
trial  profiles  were  then  constructed  and  used  in  conjunction  with  a computer 
program  for  the  direct  digital  solution  of  the  surface  integral  equations  to 
compute  the  resulting  backscattered  fields.  From  data  obtained  at  a sequence 
of  closely  spaced  frequencies,  it  was  found  possible  to  separate  the  front 
edge  contribution  from  that  of  the  rear  edge  and  traveling  waves,  and  so  choose 
the  impedance  to  minimize  each,  thereby  assuring  a broadband  performance.  Cross 
section  reductions  of  more  than  20  dB  werfe  achieved.  The  work  has  been  written 
up  as  a technical  report  (Senior,  1977)  which  will  be  sent  for  approval  of  pub- 
lication as  soon  as  the  drafting  of  the  numerous  figures  Is  complete. 
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Our  general  studies  of  the  impedance  boundary  condition  have  also  continued 
with  particular  emphasis  on  its  application  to  bodies  having  edge  and  corners. 
With  some  bodies  the  introduction  of  a non-zero  surface  impedance  has  a profound 
effect  on  the  character  of  the  solution  of  the  scattering  problem,  and  a case  in 
point  is  a plane  wave  at  skew  incidence  on  an  impedance  wedge  of  open  angle  it/2. 
If  the  impedance  of  either  surface  of  the  wedge  is  zero  or  if  the  wave  is  inci- 
dent in  a plane  perpendicular  to  the  edge,  the  solution  can  be  obtained  using 
Maliuzhinets'  (1959)  method  or,  more  simply,  by  the  method  of  images,  and 
consists  of  four  plane  waves:  the  primary  and  three  reflected  waves.  However, 
for  skew  incidence  on  a wedge  both  surfaces  of  which  have  non-zero  impedances, 
these  same  methods  are  applicable  only  if  the  (anisotropic)  impedances  satisfy 
the  compatibility  relation  of  Dybdal  et  al  (1971),  and  if  this  condition  is  not 
satisfied,  no  method  is  available  for  the  exact  solution  of  the  problem. 
Nevertheless,  it  is  possible  to  obtain  a solution  as  a Taylor  series  in  the 
impedance  n.  As  shown  by  Senior  (1978),  for  arbitrary  impedances  the  angular 
spectrum  is  no  longer  discrete,  and  the  solution  contains  an  edge  wave  whose 
amplitude  is  Ofn^)  and  whose  form  is  similar  to  that  of  the  Kottler  edge  wave 
associated  with  a 'black'  half  plane. 

In  the  area  of  low  frequency  scattering  we  have  continued  our  numerical  and 
analytical  investigations  of  the  general  polarizability  tensor  X(t).  Numerical 
values  of  the  tensor  elements  have  been  obtained  for  rectangular  dielectric 
particles  (Herrick  and  Senior,  1977)  and  the  mathematical  formulation  has  now 
been  applied  to  hexagonal  crystals  as  well.  As  in  the  case  of  rectangular 
particles,  the  tensor  is  diagonal  and  the  integral  equations  which  must  be 
solved  to  determine  the  tensor  elements  are  singular  only  along  the  edges  of 
the  crystal.  This  simplifies  the  numerical  solution  of  the  integral  equations 
by  the  moment  method  and  again  it  is  found  that  the  kernels  of  the  equations 
can  be  integrated  analytically  over  each  sampling  cell.  A computer  program  has 
been  written  to  compute  X(t)  for  a hexagonal  crystal  for  arbitrary  values  of 
T and  the  ratio  of  the  axial  length  to  the  cross  section.  However,  analysis  of 
the  numerical  results  is  not  yet  complete.  In  order  to  compute  the  tensor 
elements  for  lossy  dielectric  materials,  modified  programs  have  been  written 
which  allow  complex  values  of  t. 
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■ .1  several  aspects  of  the  mathematical 

in  addition  to  the  humenc  ^ relatiohs  which  were  assumed 

formulattoh  ha.e  beeh  examined.  The  y 

amoh,  the  Potential  functions  (see  for  examp^ 

nave  now  Poo"  P™^  \ „„  3„,,,ce  of  the  body,  are  odd 

integral  eguations.  i.e..  the  P ,„^„„„geheous  term  is  an  odd  (even) 

(even)  functions  of  the  coordinate  x,  if  the  g ^ ^ 

function  of  xi  and  the  body  the  potentials  given  by  the 

proved  that,  ^ ° ^ ,red  aero  induced  charge  condition.  In  the 

a"c  body,  such  a proof  can  easily  be  constructed  using  the 

symmetry  of  the  potentials. 
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SKEW  INCIDENCE  ON  A RIGHT-ANGLED  IMPEDANCE  WEDGE 

T.B.A.  Senior 
Radiation  Laboratory 
The  University  of  Michigan 
Ann  Arbor,  Michigan  48109 


ABSTRACT 

A problem  which  is  relevant  to  propagation  in  a rectangular  waveguide  with 
imperfectly  conducting  walls  is  a plane  wave  at  skew  incidence  on  a right-angled 
wedge  whose  surfaces  have  non-zero  impedances  either  scalar  or  tensor.  The 
solution  is  obtained  as  a power  series  in  the  impedance  n through  terms  0(n  )• 
Its  properties  are  discussed  and  it  is  shown  that  unless  an  impedance  compati- 
bility relation  is  satisfied  the  solution  contains  an  edge  wave  whose  amplitude 
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is  proportional  to  n and  whose  form  is  similar  to  the  wave  emanating  from  the 
edge  of  a 'black'  half  plane. 


INTRODUCTION 


In  the  analysis  of  scattering  problems  the  impedance  boundary  condition  is 
a valuable  tool  for  simulating  the  effects  of  surface  imperfections  in  the  form 
of  finite  conductivity  (Senior,  1959a)  or  roughness  (Senior,  1959b).  The  effects 
are  represented  via  a surface  impedance,  and  the  results  have  proved  meaningful 
even  in  the  case  of  bodies  having  corners  and  edges.  For  some  geometries  the 
scattering  problem  can  now  be  solved  using  only  a trivial  extension  of  the  method 
applicable  to  the  perfectly  conducting  (smooth)  surface,  but  with  others  the 
introduction  of  a non-zero  impedance  has  a profound  effect  on  the  analysis  and 
the  character  of  the  solution  itself.  An  example  of  the  latter  type  is  a plane 
wave  at  skew  incidence  on  an  impedance  wedge  of  open  angle  Tf/2. 

If  the  impedance  of  either  surface  of  the  wedge  is  zero  or  if  the  wave  is 
incident  in  a plane  perpendicular  to  the  edge,  the  solution  can  be  obtained 
using  Maliuzhinets*  (1959)  method  or,  more  simply,  by  the  method  of  images, 
and  consists  of  four  plane  waves:  the  primary  and  three  reflected  waves.  How- 
ever, for  skew  incidence  on  a wedge  both  surfaces  of  which  have  non-zero  im- 
pedances, these  same  methods  are  applicable  only  if  the  (anisotropic)  impedances 
satisfy  the  compatibility  relation  of  Dybdal  et  al  (1971),  and  if  this  condition 
is  not  satisfied,  no  method  is  available  for  the  exact  solution  of  the  problem. 
Nevertheless,  it  is  possible  to  obtain  a solution  as  a Taylor  series  in  the 
impedance  n and  this  is  determined  to  0(n  )•  For  arbitrary  impedances  the 
angular  spectrum  is  no  longer  discrete  and  the  solution  contains  an  edge  wave 
similar  to  that  obtained  by  Kottler  (1923)  in  his  analysis  of  diffraction  by  a 
'black'  half  plane.  The  results  are  relevant  to  propagation  in  a rectangular 
waveguide  having  two  adjacent  walls  imperfectly  conducting. 


2 


FORMULATION 

In  terms  of  the  Cartesian  coordinates  (x,  y,  z)  the  equations  of  the  two 
faces  of  the  wedge  are  y = 0,  x ^ 0 and  x = 0,  y ^ 0 (see  Figure  1).  At  each 
face  the  impedance  boundary  condition 

E - (n  • E)n  = Z n • n ^ H 

is  imposed,  where  n is  a unit  vector  outward  normal,  Z is  the  intrinsic  impedance 
of  free  space  and  ri  is  an  anisotropic  (tensor)  normalized  impedance.  On  y ■ 0 

n = XX  + n2  zz 

and  the  boundary  conditions  are  therefore 

E,  -n,  ZH^.  E^'-njZH,  . (1) 

whereas  on  x = 0 

n = 03  yy  + n4 

giving 

Ey  = - na  z Hz  . E^  = n4  z Hy  . (2) 

The  incident  field  is  a plane  wave  of  arbitrary  polarization  incident  in 
a direction  making  an  angle  b with  the  negative  z axis,  with  propagation  vector 

= - k (x  cos  a cos  B + y sin  a cos  b + z sin  b) 
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where  a time  factor  e has  been  assumed.  Since  the  surfaces  are  independent 
of  z,  the  total  field  will  depend  on  z only  through  the  factor  e , and 

can  be  expressed  in  terms  of  the  single  component  Hertz  vectors 


n = S U (X,  y)  n*  . S V (x.  y) 


as 


iU 

ax 


ikf  , 

Z H = 

- ih 

av 

ay  ’ 

X 

ax 

ik^  , 

Z H = 

- ih 

ax  ’ 

y 

ay 

au 

ay 


= x"  u , 


U and  V satisfy  the  scalar  wave  equation 

f .2  2 , 

-3-  + ^ + 

,2,2^ 

ax  ay 


Z U 


H=  0 


and 


implying 


h = k sin  6,  X = k cos  e 


h2  . ^ k^  . 


The  boundary  conditions  on  U and  V are  then 

h a . . fk  ^ * -L  ,2  ] u 
ax  [ 3y  09  j 


. au  ,3^4  2 

h— = k— +in,  X 
ax  ay  1 


V . 


n T~  * 

ay 


k |-  + in,  x^l  V . h = 

ax  3 I * ay 


ax  n4 


(3) 


(4) 

(5) 


on  y » 0 and  x =»  0 respectively,  and  serve  to  couple  U and  V unless  h = 0 
(incidence  in  a plane  perpendicular  to  the  edge),  n2  “ 04  * 0 or  r>^  “ 03  ' 


4 


To  specify  the  incident  field  we  choose 


ui  = A = B ^ (6) 

where  a factor  ^ has  been  suppressed.  A and  B are  arbitrary  con- 

stants and  c = cos  a,  s = sin  a.  Since  a plane  wave  cannot  propagate  parallel 
to  the  surfaces,  it  is  necessary  to  assume  that  c,  s > 0 implying  0 < a < tt/2, 
and  X f 0.  The  task  is  to  find  the  total  (incident  plus  scattered)  fields  in 
accordance  with  (6)  and  the  boundary  conditions  (4)  and  (5). 


PLANE  WAVE  SOLUTION 


We  first  seek  a plane  wave  solution  in  the  form 

U = A + A^  + A^  + A3 

(7) 


V = B + B^  + B2  + 83 


(8) 

as  suggested  by  the  image  method.  The  last  three  terms  in  each  of  (7)  and  (8) 


constitute  the  scattered  field,  and  the  unknown  amplitudes  A^  and  B^,  i = 1,  2,  3, 
must  be  chosen  to  satisfy  the  boundary  conditions  (4)  and  (5).  From  the  con- 
ditions on  y = 0 

he  (A  + A^)  = (ks  - Xn^ ) B - (ks  + Xn^)  B^  (9) 

he  (A2  + A3)  = - (ks  - xn^)  B2  + (ks  + Xn^)  B3  (10) 

he  (B  + B^)  = - (ks  - X/n2)  A + (ks  + X/ri2)  (11) 

he  (B2  + B3)  = (ks  - x/n2)  ^2  ' ^'^'’2^  ^3  ' 
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Similarly,  from  the  conditions  on  x = 0 


hs  (A  + A^)  = - (kc  - Xn^)  B + (kc  + An^)  (13) 

(A^  + A^)  = (kc  - An^)  B^  - (kc  + An^)  B^  (14) 

hs  (B  + 82)  = (kc  - A/n^)  A - (kc  + A/n^)  (15) 

hs  (B^  + B3)  = - (kc  - A/n^)  A^  + (kc  + A/n^)  A3  (16) 


yielding  a total  of  8 equations  for  the  6 unknowns  A^  and  B^.  From  (9)  and 


(11)  A-j  and  can  be  expressed  in  terms  of  A and  B: 

rA^  = - {r  - 2ks  (ks  + An^)}  A + 2 hksc  B (17) 

rB^  = - 2 hksc  A - {r-2ks  (ks  + A/n2)>  B (18) 

where 

r = h^  c^  + (ks  + Ani)(ks  + \/^^)  . (19) 

In  addition,  from  (13)  and  (15) 

r'  A2  = - {?'  - 2kc  (kc  + An3)}  A - 2 hksc  B (20) 

r'  B2  = 2 hksc  A - {r*  - 2kc  (kc  + A/n4)}  B (21) 

where 

r'  = h^  s^  + (kc  + An3)(kc  + A/n4)  (22) 


completing  the  specification  of  all  the  amplitudes  save  A3  and  B3. 

From  (10)  and  (12),  however,  A3  and  B3  can  be  expressed  in  terms  of 
^2  3nd  B2,  and  using  then  (20)  and  (21),  we  have 


rr'  A3  » Qr  - 2ks  (ks  + An^)}{r'  - 2kc  (kc  + An3)}  - 4 hVsV^A 

+ 2 hksc[]r  - 2ks  (ks  + An^)  + r*  - 2kc  (kc  + A/n^)  Jb  (23) 

6 il 


rr'  = - 2 hksc  Qr  - 2ks  (ks  + A/n2)  + r'  - 2kc  (kc  + Xn2)^A 

+ [{r  - 2ks  (ks  + \/T)^)}{r'  - 2kc  (kc  + x/^^))  - 4 h^k^s^c^Js. 

(24) 

Also,  from  (14)  and  (16)  with  (17)  and  (18) 


rr 


' Aj  = Qr  - 2ks  (ks  + Xn^)}{r'  - 2kc  (kc  + Xn3)}  - 4 h^k^s^c^^A 

- 2 hksc|^r  - 2ks  (ks  + X/n2)  + r'  - 2kc  (kc  + Xn3)  Jb  (25) 


rr'  B3  = 2 hksc  j^r  - 2ks  (ks  + Xn-j)  + r'  - 2kc  (kc  + x/n^)^A 

+ [{r  - 2ks  (ks  + x/n2)}{r'  - 2kc  (kc  + x/n^))  - 4 h^k^s^c^i 

(26) 

and  the  results  are  not  necessarily  equivalent.  Taking  first  the  expressions 
(23)  and  (25)  for  rr'  A^,  the  coefficients  of  A are  identical  but  those  of 
B can  be  written  as 

- 2 hksc  [[kx  {s  (n^  - l/n2)  “ c (n3  - 1/04)}  TX^(^^  + ^ - l]ll 

(27) 

with  the  upper  and  lower  signs  for  (23)  and  (25)  respectively.  Since 

X,  s,  c ^ 0,  the  coefficients  are  the  same  if  and  only  if  h = 0 corresponding 

to  incidence  in  a plane  perpendicular  to  the  edge,  or 

Hi  n-j 

_L  + _i  . 1 = 0 . (28) 

n2  04 

Similarly,  in  the  expressions  (24)  and  (26)  for  rr'  83,  the  coefficients  of 
B are  identical  but  those  of  A are  now  given  by  (27),  leading  again  to  the 
requirement  that  (28)  be  satisfied  to  make  them  the  same.  Equation  (28)  is 
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the  impedance  compatibility  relation  obtained  by  Dybdal  et  al  (1971)  as  the 
condition  for  the  existence  of  discrete  modes  in  a rectangular  waveguide  whose 
opposite  walls  have  the  same  anisotropic  impedances.  It  is  equivalent  to 

E • H = 0 

at  the  edge,  and  if  it  is  satisfied 

rr’  = \jT  - 2ks  (ks  + Xn^)}{r'  - 2kc  (kc  + Xn3)}  - 4 h^k^s^c^^l* 

- 2hk^  X sc  {s  (n^  - 1/02)  - c (n3  - l/n^)}  B (29) 

rr'  B3  = - 2hk^  x sc  {s  (n^  - 1/02)  - c (03  - Vn^)}  A 

+ |^{r  - 2ks  (ks  + x/n2)}{r'  - 2kc  (kc  + x/n^)}  - 4 h^k^s^c^^B. 

(30) 

The  expressions  (7)  and  (8)  for  U and  V then  constitute  a valid  solution  with 
the  amplitudes  A.  and  B.,  i = 1,  2,  3,  defined  by  (17),  (18),  (20),  (21),  (29) 
and  (30).  Since  U and  V and,  hence,  and  are  non-zero  at  the  edge,  the 
components  of  the  electric  and  magnetic  currents  perpendicular  to  the  edge  are 
continuous  from  one  face  to  the  other,  but  as  a result  of  the  boundary  con- 
ditions the  component  of  the  electric  current  parallel  to  the  edge  is  discon- 

tinuous unless  02  = and  the  same  component  of  the  magnetic  current  is 
discontinuous  unless  ni  = n^. 

Our  final  remarks  are  for  the  case  when  the  impedances  are  small.  If 
= e^n  i = 1,  2,  3,  4 
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where  the  e.  are  of  order  unity  and  n <<  1 » the  discrepancies  between  (23) 
and  (25)  and  (24)  and  (26)  are  O(n^).  It  follows  that  to  0(n)  our  plane  wave 
solution  is  valid  regardless  of  the  compatibility  relation  (28),  and  it  is 
now  of  interest  to  consider  the  form  taken  by  the  next  higher  order  term  for 
small  but  arbitrary  impedances. 

SMALL  IMPEDANCE  SOLUTION 


To  obtain  a solution  which  holds  regardless  of  the  relation  (28),  we 
assume  the  following  series  expansions  for  U and  V: 


U (x,  y)  = I o'"  Ujjj  (x,  y)  , 


V (x,  y)  = I n*"  V (x,  y) 


(31) 


m=o 


m=o 


On  substituting  into  the  boundary  conditions  (4)  and  (5)  and  equating  co 
efficients  of  like  powers  of  n,  we  find 

3V 


" ^ ^ ^ ''.-1  • ' - '2 


, m-  I . . m-  I 
k + h 


3y 


(y  = 0) 


3x 


L,  "1  . t m 

^ 3y  ^ ^ 3x 


i e,  , , i 

3 m-1  m 


- e. 


3U 


m-1 


3V, 


3x 


- h 


m-1 


3y 


(x  = 0) 


for  m ^ 0 with  U_i , V_^ 


u„  (X,  ,) 


= 0,  and  these  can  be  expressed  as 
f^l^(x)  y = 0,  x>0 


f(2)(y)  X = 0,  y > 0 


(32) 
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where  n is  the  outward  normal  and 


9*l’(x) 

g^^hy) 


y = 0,  X 1 0 
X = 0,  y 1 0 


(33) 


m ^2 


. Vl  ^ . '^m-1 

k — r—  + h 


3y 


3X 


y = 0 


(34) 


f(2)(y)  = ^ 

\ 


3U„  1 3V„  1 

k - h 


3X 


3y 


X = 0 


(35) 


g^^^(x)  = - 
^ m ' ' 


kX‘ 


3^  V 

e,  V„  , - eo  h^  - Co  kh 


3^  U, 


m-1 


-1  ''  'm-1  ’"2 


2 2 
3X  3X3y  J 


C36) 


y = 0 


g(2)(y)  = . i 


P 3^  V 1 
2 m-1 


3^  U, 


Go  X’'  V„  , - e,  h"  ^ - G.  kh 


m-1 


'3  " 'm-l 


2 '"A 

3X  3X3y  j 


. (37) 


X = 0 


The  task  has  now  been  reduced  to  the  solution  of  separate  boundary  value 

problems  for  U_  and  V . 

mm 


Given  a scalar  field  !j)(x,  y)  satisfying  the  wave  equation  (3),  the  field 
at  any  point  inside  and  on  a closed  contour  L can  be  written  as 


<(-(x.  y)  = [-f  - G ^ ds* 


(38) 


where  G is  a two-dimensional  Green's  function.  In  the  case  iji  = it  is 
convenient  to  choose  G = G^: 
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with 


V/2 

*^1,2" 

x*)^  + 

(y 

; y')^ 

] 

r 

^1/2 

*^3,  4 " 

|^(x  + 

x')'^  + 

(y 

+ y')^ 

) 

is  the 

Hankel 

function 

of  the 

first 

G vanishes  on  the  surfaces  x = 0 and  y = 0.  On  applying  (38)  to  the  region 


bounded  by  the  lines  y = 0 (x  ^ 0) , x = 0 (y  ^ 0)  and  a quarter  circle  at 
infinity,  and  imposing  the  boundary  conditions  (34)  and  (35),  we  have 

5 

3y' 


dx'  + 


oi 

/ 


m ' 9x' 


y‘=0 


dy- 


x‘=0 


since  the  quarter  circle  yields  no  contribution.  Hence 


U 


m 


(x,  y)  = 


00 

r '1 

J m ' ] 9y  0 

0 ^ 

CSJ 

+ 

CNJ 

X 

1 

X 

(1) 

0 


X + x' ) 


_3 

3X 


,(1) 


/(y  + y')^  + x^ 


dy 


(39) 


from  which  U can  be  computed  knowing  f^^^(x)  and  f^^^(y). 
m ^ m m 

For  the  partial  field  the  appropriate  choice  of  6 is  G = where 

and  using  this  and  proceeding  in  an  analogous  manner  to  the  above,  we  find 


n 


which  specifies  in  terms  of  and  g^^^(y). 

The  incident  plane  wave  is  given  by  (6),  and  since  the  boundary  conditions 
on  the  zero  order  fields  are 


on  y = 0,  X ^ 0 
on  X = 0,  y ^ 0 

it  follows  immediately  that 

U (x.  y)  = A(e"''^^^^  ‘ ^^0 

(41) 

Vq  (x,  y)  = B ^ ^"0 

(42) 

We  can  now  determine  the  boundary  values  for  the  first  order  fields.  From 
(34)  and  (35)  wi th  m = 1 , 

f^l^(x)  = - (ks  A + he  B)  sin  x xc  , 
f'^  (y)  “ — (kc  A - hs  B)  sin  xys  , 


!!c 

ay 

!!c 

ax 


= 0 


= 0 
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and  when  substituted  into  (39),  we  have  (see  A. 3) 


l),  (X,  y)  = - ^ (ks  A ♦ he 

. (kc  A - hs  B)  . e *»<«  - . (43) 

Similarly,  for  , 

g^l^(x)  = - ^^2  ^ ■*■  ^^1  ■*■  ^2  ’ 

g^^^(y)  = ^1^4  khsc  A - (e^  + £4  h^s^)  bJ  cos  xys  , 

and  when  substituted  into  (40),  the  result  is  (see  A, 5) 

V,  (X.  y)  > - L khsc  A t (c, 

* kfe  (‘4  » - <=3  M • 

(44) 


These  serve  to  specify  the  boundary  values  for  the  second  order  fields. 
From  (34)  and  (35)  with  m = 2, 


f^l^(x)  = 


4iee 


kx‘s 


(■ 


ke,s  (k^  - x^c^)  A + he  [e^  X^  + g2  - X^c^)]  B )sin  xxc 


0 


+ ^ 1^:4  ksc  A + h (e3  - £48^)  bJ 


kc.c  (k^  - x^s^)  A - hs  [e,  X^  + £.(k^  - x^s^)]  B )sin  xys 


^ kx'^c  L 

+ -j^  |g2  ksc  A - h (ei  - £2^^)  6 


3 A £44 

ixys 
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giving 


U2  (x,  y)  = 


2e2 

kx^s 


^k£2S  (k^  - x^c^)  A + he  [e^  X^  + (k^  - x^c^)] 


_ p-i^(xc  - ys)}+  — 2-  p-4 
k X c 


)J+  — ^ ^kc^c  (k^  - x^s^)  A - hs  [£3  x^  + £^  (k^  - x^s^)] 
^ei^(xc  + ys)  . 3ix(xc  - ys)j  , ksc  A + h (£3  - £^52)  ej  (x,  y;  a) 

+ -j^  ^£2  ksc  A - h (£^  - e2C^)  B^  (y,  x;  - a) 

where  the  function  is  defined  in  the  Appendix.  As  shown  there,  however, 

(y.  x;  j - a)  = ^ (x,  y;  a)  , 


and  U2  can  therefore  be  written  as 


U2  (x,  y) 


‘■^2 

kx^s 


^ke^s  (k^  - x^c^)  A + he  [e^  x^  + E2(k^  - X^c^)] 

- ^ *-^3  ^ 

|e^A(xc  + ys)  . ^ix(xc  - ys)j^  ksc  A - h (e^  - £20^)  sj 


4h 


+ ^ ■*■  ^2^3  ' ^2*^4^  ® ‘^1 


(45) 
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Similarly 


g^2^(x)  = 


^£2  + ^2  ^+[(^1  + ^2  * ('^2 


k"xs 


cos  X xc 


ix  xc 


g^2^{y) 


4i 

Sc 


J-  khsc  [£3  x^  + £^  (k^  - x^s^)]  A - [(£3  X^  + E^  h^s^)^  <■  (e^  khsc)^]  sj 


• cos  X ys 


+ ^^|e2  khsc  (£3  - e^s^)A  + [(£3  - e^s^)(£^  x^  + e^  h^c^)  + e^e^  k^s^]  bJ 


• e’'  . 

On  substituting  these  into  (40)  and  using  (A. 4),  (A. 12),  (A. 14)  and  (A. 15), 
we  find 


V2  (x,  y)  = 

-g- 1 2 ^£2  * ^2  ^ * ^2  + (^2  khsc)^] 

. ^ei^(xc  + ys)  ^ g-ix(xc  - ys)^ 
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waveguide  having  two  adjacent  walls  Imperfectly  conducting. 

2 


khsc  [€3  > E4  (k^  - x^s^)]  A - [{€3  + (e^  khsc)^]  E 

A C Vh  .. 

. |^eiA(xc  + ys)  ^ g1x(xc  - ys)J 

^£3  khsc  (£3  - e^s^)  A + [(£3  - £^s^)(£^  x^  + £3  h^c^)  + £^£^  k^s^] 


k‘-sc 


4h 


1x{xc  + ys) 


” k ^^1^4  ^^2^3  "^2^4^  ^ *^2 
with  J2  defined  in  the  Appendix. 


(46) 


Since  A = B = 0 implies  a null  solution,  the  last  term  in  each  of  (45)  and 
(46)  vanishes  only  if  the  impedance  compatibility  relation  (28)  is  satisfied. 

U2  and  V2  then  consist  of  three  plane  wave  contributions,  and  the  resulting 
expressions  for 

U = Ug  + n U]  + U2  . V = Vjj  + n V2 


O 

are  identical  to  the  expansions  through  terms  0(n  ) of  the  solution  derived  in 
the  previous  section.  If  desired,  higher  order  terms  can  also  be  obtained,  and 
involve  only  the  plane  waves  already  present  in  U2  and  V2.  If  (28)  is  not 
satisfied,  however,  the  character  of  the  solution  changes  considerably.  U2 
and  V2  now  contain  terms  proportional  to  and  J2  respectively,  and  it  is  of 
interest  to  examine  the  solution  in  this  more  general  case. 


DISCUSSION 

The  functions  (x,  y;  o)  and  J2  (x,  y;  o)  satisfy  respectively  the 
integral  equations 
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f 


3G. 


f (x.  y)  = / f (x',  y')  ^ds' 


and 


g (x,  y)  = 


■/ 


3 

'h  3n' 


Gu  g (x’ , y' ) ds' 


where  l is  the  contour  consisting  of  the  straight  line  segments  y = 0,  x ^ 0 
and  X = 0,  y ^ 0.  Each  can  be  expressed  (see  Appendix)  in  terms  of  the  function 


( 

Mx  . 4-)  = j" 


sin  !p ircosh  t 

cosh  t + cos  ifi 


(47) 


which  occurs  in  Kottler's  (1923)  solution  for  the  diffraction  of  a plane  wave 
by  a perfectly  'black'  half  plane.  The  integral  is  absolutely  convergent  if 
Im.  T > 0 and  -ir  < 4;  < w,  and  has  been  discussed  at  length  by  Watson  (1938), 
Copson  and  Ferrar  (1938)  and  Erdelyi  (1938a,  b).  It  is  similar  to  the  Fresnel 
integral  and  can  be  expanded  in  a series  of  'cut'  Hankel  functions  (Copson  and 
Ferrar,  1938).  In  addition, 

F (,  , ♦)  » [sin  (48) 

0 

(Erdelyi,  1938a),  constituting  a finite  integral  representation,  but  no  analytic 
evaluation  of  the  integral  has  yet  been  achieved.  For  t » 1 the  stationary 
phase  approximation  to  (47)  is, 

and  for  t « 1 (48)  shows 


F (x  , 4»)  “ ^ + 0 (t  log  t). 
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k. 

From  the  expressions  (A. 9)  and  (A. 13)  for  J-j  and  respectively  in  terms 
of  F(t  , III),  it  follows  immediately  that  on  the  line  segment  y = 0,  x > 0 


1 - -U  X cos  a 
Jl  - e 


aJ 

3y 


2 _ ..  ^ix  X cos 

1 X sin  a e 


whereas  on  the  segment  x = 0,  y > 0 


and  are  also  closely  related  and  if  (x,  y;  a)  = K (e  , «)  where  e is 
the  cylindrical  polar  coordinate,  (x,  y;  a)  = K (ir/a  - a,  Tr/2  - e).  For 
Xp  » 1 and  9 bounded  away  from  a 


(x,  y;  a) 


giXp  cos  (a 


sin  29 

cos  2d  - cos  2a 


^2  ’ y’’ 


giXp  cos  (a  - 9)  H 


i (xp  + Tr/4) 


sin  29 , 

cos  20  - cos  2a 


where  H is  the  unit  step  function,  showing  that  at  large  distances  from  the 
edge,  both  functions  have  the  appearance  of  an  edge  wave.  For  Xp  « 1 , however. 


(x,  y;  a)  ~ 1 - 29/Tr  , (x,  y;  a)  ~ 2a/T(  . 

The  partial  field  V2  is  therefore  continuous  from  one  surface  of  the  wedge  to 
the  other,  but  has  a finite  jump  discontinuity. 

Partial  fields  and  of  higher  order  than  the  second  can  be  obtained 
by  the  same  method  used  to  determine  II2  and  V2,  and  these  turn  out  to  be  singular 
at  the  edge,  with  the  order  of  the  singularity  Increasing  with  m.  They  also 

i 
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involve  functions  which,  though  similar  to  J-j  and  become  more  complicated 
as  m increases,  and  for  this  reason  it  does  not  seem  likely  that  the  field 
expansions  can  be  summed  to  arrive  at  an  exact  closed  form  solution  of  the 
wedge  problem.  Nevertheless,  it  is  evident  that  like  the  partial  fields 
and  for  m ^ 2,  the  exact  solution  must  have  a continuous  angular  spectrum. 
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APPENDIX:  SOME  INTEGRAL  EVALUATIONS 


There  are  two  types  of  integral  that  arise  in  connection  with  the  partial 
fields  U-j  and  V^ . The  first  is 

= Jsin  (A  x'  c)  H^^^a/(x  - x')^  +7)  - ^ H^^^a/{x  + x')^  + y^) 


(A.l) 


which  can  be  written  as 


00 

= J sin  (A  X'  c)  3^H^J^(a/(x  - x')2  +y2) 


We  now  insert  the  Fourier  integral  representation 


, ^2  _ 1 (x  - x')  + iy  '/k  - u)  doi 


H'q^(a/(x  - x')^  + y^)  = ^ e 


- J- 


(A. 2) 


valid  for  y ^ 0,  where  A is  assumed  to  have  a small  positive  imaginary  part. 
On  reversing  the  order  of  integration  and  using  the  fact  that 


L-i  (w  T Ac)  X’  = 2Tr  6 {<.  T Ac) 


where  6 is  the  Dirac  delta  function,  we  have 


The  second  integral  is 


(A. 3) 


I2  = AS 


J"cos  (a  x'  c)  ^H^q^(a/(x  - x' )^  + y^)  + H^q^(a/(x  + x')^  + y^^dx' 


(A. 4) 
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and  in  a similar  manner  it  can  be  shown  that 


I = (e’^  + e'"' 


(A. 5) 


For  integrals  that  differ  from  the  above  in  having  only  single  exponentials 
in  the  integrands  the  evaluation  is  more  difficult.  Consider 


(x»  y»  ®)  “ 

(A. 6) 

where  0 < a < tt.  If  we  again  insert  the  Fourier  integral  representation  (A. 2), 
the  x'  integration  can  be  carried  out  immediately  to  give 


gi  (x  cos  a + (o)  x'  ^ 


X cos  C(  T u) 


and  hence 


J]  (x 


2 2 2 
0)  - X COS  a 


iwx  ^ iy 


We  now  make  the  substitution  w = x cos  8 and  write 


X = p cos  0 , 


with  0 ^ 0 £ 7t/2,  so  that 


y = p sin  0 


J,  (X,  y;  a)-  -if  - ^ dS 

” J cos^  6 - cos'^  “ 

C 

where  C is  the  path  shown  in  Figure  A-1.  The  stationary  phase  point  is  8 = e , 
and  in  a deformation  of  the  path  into  the  stationary  phase  contour  $(0)  the  pole 


at  8 = a is  captured  if  e < a.  Hence 


J,  (x,  y;  a)  = e’^"^  H (a  - 9) 


_ i / sin  (8  •*•  e 

”s(0)  <«*  ®> 


) cos  ( 8 -I-  9) 


COS^  a 


giXp  COsX^jj 


(A. 7) 


where  H denotes  the  unit  step  function. 

The  path  S(0)  is  symmetrical  about  the  origin  of  the  X plane,  and  the  method 
employed  by  Clemmow  (1951,  p.  293)  can  now  be  used  to  reduce  the  integral  to  a 
more  elementary  form.  Replacing  X by  -X  and  then  adding  the  resulting  integral 
to  that  in  (A. 7),  we  find  after  some  manipulation 

(x,  y;  a)  = e'^'’  ‘ H (a  - 9)  + 

sin  (9  g)  sin  (9  - a)  sin  (9  - a) 

cos  i + cos  (9  + a)  ” cos  X + COS  (9  - a)  “*■  COS  5 - COS  (6  - a) 

• e’^*^  dX  . (A. 8) 


sin  (9  + g) 

COS  X - COS  (9  + a) 


Each  cf  these  integrals  is  only  a special  case  of  the  Kottler  function  F (t,  ij;) 
defined  in  (47).  A simple  change  in  the  variable  of  integration  shows 


F (Xp  , ii)) 


4ir  I c 

do) 


sin  It 


,Up  cosX 


cos  X + COS  ill 


dX 


and  with  the  restriction  - ir  < iji  < tf,  the  first  three  integrals  in  (A. 8)  cor- 
respond to  the  identifications  \j/  = TT-9-a,  9 + a and  e - a respectively. 
For  the  fourth  integral  i(/  = tit-e  + a according  as  9 < a respectively,  and 
since  (48)  implies 


F (Xp  , ff  - 9 + g)  - F(Xp,  - IT  - 0 + g)  ® 
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it  follows  that 


(x,  y;  a)  = 


F (xp  , TT  - 9 - a)  - F (\p  , 6 + a)  - F (Xp  , 0 - a)  + F (Xp  , tt  - 0 + a). 

(A. 9) 

On  interchanging  x and  y in  (A. 6)  and,  at  the  same  time,  replacing  a by 
Tt/2  - a,  we  have 


(y.  x;  tt/2  - a)  = 

00 

- J H^J^hx/(y  - y')^  + x^)  - H^J^x/(y  + y')^  + x^)j  dy* 

O 

= F (Xp  , 0 + a)  " F (Xp  , TT  - 0 - a)  - F (Xp  , a - 0)  + F (Xp  , tt  + 0 - a) 

(A. 10) 


and  hence,  from  (48), 


J]  (x,  y;  a)  + (y,  x;  tt/2  - a)  = e^^^  (a  - 0)  ^ 


(A. 11) 


The  only  integral  still  to  be  evaluated  is  the  analogue  of  (A. 4)  with  the 
cosine  factor  replaced  by  a single  exponential.  If 


J2  y>  “)  “ 


LAin  a J* gixx'  cosa|J,(l ) - x')^  + y^)  + H^g^x/(x  + x')^  + y^)J  dx'. 


(A. 12) 
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insertion  of  the  Fourier  integral  representation  of  the  Hankel  function  gives 
J2  y>  a)  = 


/•  ^ 

ix^  / sin  m 

■ ~J^ 

•/  u)  - A 


cos  a luX  + ly 


+ iy/x^  - 


COS^  a 


2 

VA  • (jO 


IXp 


H (a  - e)  - i r 

' ■’s(O) 


sin  g cos  a 

2 2 
cos  {X  + 9)  - cos  a 


,ixp  cosj 


dX  , 


and  by  a process  similar  to  that  used  above,  we  find 


Jo  • y»  o)  ~ 


gixp  cos  (a -9)  H („  - e)  - 


-i/  (, 


sin  (9  g)  sin  (e  g) 

cos  X - cos  (0  + a)  " cos  ^ + cos  (0  + a) 


+ sin  (9  - g) sin  (e  - a)  . cos 3 

cos  6 + cos  (9  - a)  “ CoTX  - cos  (9  - a)  / ® 


Hence 


(x,  y;  a)  = 


- F (Xp,  TT  - 0 - a)  + F (Xp,  0 + a)  - F (Xp,  0 - a)  + F (Xp,  tt  - 0 + a). 

(A. 13) 


With  9 and  a replaced  by  tt/2  - 9 and  it/2  - a respectively, 

^2  (y*  - a)  = 

Lc^f  sin  i ,■ )?  * 

o ''  > 

= - F (xp,  0 + a)  + F (Xp,  It  - 0 - a)  - F (Xp,  o - 0)  + F (Xp,  tr  + 0 - a)  , 
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(A. 14) 


FIGURE  CAPTIONS 


Figure  1:  The  geometry. 

Figure  A-1:  Path  of  integration  in  complex  a plane. 


